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Abstract 



The scattering of two off-shell photons is an infrared-safe process in QCD. For 
photon virtualities Q in the range of a few GeV, accessible at LEPII, power- 
behaved contributions in (Aqcd/Q)™ to the total cross section may become 
non-negligible. Based on a dispersion relation for the running coupling, we 
discuss these contributions and calculate the coefficients of the leading power 
correction for transversely polarized and longitudinally polarized virtual pho- 
tons. 



1 



Typeset using REVTeX 



Two-photon collisions provide one of the dominant modes in experiments at high-energy 
e ± e~ colliders. If the photons are sufficiently far off shell, the process is dominated by 
short-distance QCD interactions. A survey of the QCD studies in photon-photon scattering 
that are currently being carried out at LEPII may be found in Ref. 0. In particular, the 
total cross section for scattering two off-shell photons is an infrared-safe observable. The 
role of this property has recently been emphasized in the context of investigations of the 
high-energy limit of QCD [J2J. Predictions for (7(7*7*) can be computed in perturbation 
theory and tested as a function of the photon virtualities. For virtualities Q of the order of 
a few GeV, contributions to o suppressed by powers of Aq CD /Q may become non- negligible. 
These contributions are the subject of this paper. 

A systematic approach to the calculation of power-behaved corrections to hard processes 
is an open problem in QCD. For cases in which an operator product expansion is applicable, 
this provides a general framework to classify higher-twist contributions. However, so far this 
has proved to be of only limited practical use. On the other hand, there have been efforts to 
develop methods that allow one to derive estimates of power corrections to a given observable 
from the study of the infrared behavior of its perturbative series [||] . In this context different 
techniques have been proposed over the past few years and applied to a variety of hard 
processes (for recent reviews see, for instance, Ref. |3J). The basic observation underlying 
these methods is that the factorial growth of the coefficients of the QCD perturbation series 
in large orders gives rise to ambiguities in the perturbative predictions which are proportional 
to power-behaved contributions. These ambiguities can be interpreted as being due to an 
artificial separation between short-distance and long-distance physics in the perturbative 
treatment. From the requirement that they must cancel in the physical cross sections once 
higher orders as well as nonperturbative contributions are included, one is able to derive 
information on the structure of the power correction. 

The physical origin of these power-like contributions is an infrared one. They are as- 
sociated with the loop integrations over the regions of small momenta in Feynman graphs. 
Based on this, the authors of Refs. [|]|| have proposed a dispersion relation for the QCD 
running coupling in order to relate power-like corrections to the behavior of the coupling 
at small momentum scales. In this paper, we will use this dispersion relation to analyze 
the exchange of gluons in high-energy photon-photon interactions. This will enable us to 
identify the leading power correction to the total photon-photon cross section. 

We will discuss the scattering of two off-shell (spacelike) photons with momenta qA and 
Qb, 7*(<?a) + 1*(<1b) hadrons, with the virtualities q\ = —Q\ and q 2 B = —Q 2 B being 
large compared to Aq CD . We will focus on the region where the center-of-mass energy 

y/s = yj (qA + qB) 2 is much larger than Qa and Q B . In this region questions related to power- 
like corrections become especially important, because power-like corrections are expected to 
be associated with the mechanism that unitarizes the total cross section at asymptotic 
energies. We will thus start with the large-s form of the perturbative cross section, in 
which terms that fall like Q 2 /s are neglected. In the Born approximation, the corresponding 
amplitude is given by graphs with exchange of one gluon between two quark-antiquark pairs 
created by the virtual photons (Fig. [I]). In this approximation the total cross section has 
the structure [Q] 
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G(k 2 ,Q 2 A )G(k 2 ,Q 2 B ) 



(1) 



Here d 2 k denotes the integration over the transverse momentum flowing in the gluon line. 

The factors 1/ (k 2 ^j come from the gluon propagators. The factors G are each proportional 

to a s , 



G = a s g(k 2 /Q 2 



(2) 



and describe the coupling of the gluon to the q q system. The explicit form of the functions 
g(k 2 /Q 2 ) depends on the photon polarization. In what follows we will discuss first the cross 
section averaged over the two transverse polarizations and then we will extend our results 
to the case of the longitudinal polarization. 
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FIG. 1. One of the graphs contributing to the Born amplitude for the process 7* +7* 
at large s. 
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The cross section o"o has a scaling behavior in the photon virtualities of the type 
a oc 1/Q 2 , modulated by logarithmic scaling-violation factors in the ratio of the two virtu- 
alities. To evaluate the corrections to this behavior that are suppressed by powers of the 
photon virtualities, we begin by introducing the running coupling a s (k 2 ) for each one of the 
factors G in the integral (|l[). Following the method of Refs. f|||, we consider the dispersive 
representation of the QCD coupling in terms of the spectral density p s , 



a s (k 2 ) = - [ 
Jo 



00 djjL 2 



/o y 2 + k 2 

and we define the effective coupling a e fr according to the relation 
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We thus have 



a,(k 2 ) = k 2 f 
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The effective coupling a e ff differs from the perturbative coupling in the infrared region. The 
form of a c g is determined by nonperturbative physics. As we shall see, the main point 
of this approach is that, if this form can be assumed to be universal, then power-behaved 
contributions to the cross section can be parametrized in terms of moments of the effective 
coupling variation at low transverse momentum scales. 

Replacing the strong coupling in Eq. flU) by Eq. @ yields 

Z7T JO fi A JO fi B 

with the function being given by 

poo A Xp' 

4>{<? A ,<?B\A,&)=A& I {k 2 ~ 2 )2 (k 2 + 2 )2 9^ IQb) • (7) 

Eq. (H) is written in terms of two dispersion variables (jla, A*b, corresponding to the fact that 
our starting amplitude (Fig. [I]) is a two-loop contribution in QCD. Introducing the variation 
5a e f[ of the strong coupling at low scales, we decompose a e s as follows 

a e ff(/i 2 ) = «PT(/i 2 ) + £a eff (/i 2 ) , (8) 

where « P t denotes the form of the coupling in perturbation theory. The decomposition (H) 
gives rise to terms in a PT x a PT , a PT x 5a e ff, and £a e ff x fi&efi in Eq- @- 

The integration of the terms in a P x x a P x is dominated by values of the dispersion 
variables of the order of the hard scales, Qa and Qb, as we will see below. Thus, these 
terms give rise to the standard leading result for the cross section in perturbation theory. 
The terms involving £a e fr, on the other hand, probe the behavior of the function at small 
values of the dispersion variables, because the distribution of 5a e g is concentrated at small 
scales. These terms are responsible for power-behaved corrections to a. We will see that in 
the case of the leading power correction the quadratic term in 6a e s does not contribute and 
the correction comes entirely from the linear term in <5a e fr- 

To study the form of the function <fi at small /i 2 , we use the explicit expression for g. For 
transversely polarized photons, this is given by 

9{kW) = 8 a (E eij f „ , P« / „ A P(A) - (j _ - (i - - - , (9) 

where a is the electromagnetic fine structure, e q is the quark electric charge in units of 
e = v47ra, and P is the vector — > fermion splitting function 

Hx) = *±<£*r . (10) 

From Eqs. (0) and ([|) we may observe that for small fi the function vanishes. This is 
associated with the fact that the 7* 7* cross section is an infrared safe quantity in pertur- 
bation theory (see Eq. (jl|)), that is, the factors g vanish like k 2 /Q 2 (times logarithms) as 
k 2 — >• 0. For large \± the function <fi also vanishes. This is associated with the fact that the 
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FIG. 2. The dependence of the function <p on one of the dispersion variables (^b) for dif- 
ferent values of the other one (ha)- We take Qa = Qb and plot the rescaled function 
$ = QaQb 4>/ (8 a £,e*) 2 versus h%/(QaQb)- 



graph in Fig. [I] is well-behaved in the ultraviolet, that is, the factors g do not grow more 
than logarithmically at large k 2 . 

We now evaluate the function <fr numerically. In Fig. |2] we plot the /i^-dependence of 
the dimensionless function QaQb4> f° r fixed (small) values of [Aa, working at Qa = Qb- 
First, we note that <fi peaks in the region where the ratio /za/Vb of the dispersion variables 
is of order 1 and their product is of order QaQb- Then the leading contribution to the 
term apx x ctpx from the integral (H) is obtained by pulling out of the integral the factors 
of ctpT evaluated at a scale of order Qa Qb- We thus recover the leading perturbative result 
(HI). Second, we consider the term <5a eff (ji A ) a PT (/i 2 3 ). We see from Fig. || that for small 

peaks at a value ~p B proportional to \JQaQb and to 



[LA the function peaks at a value ~pr B proportional to \[QaQb and to \j[i A - Then we 
approximately perform the /i^-integration by pulling out of the integral the factor of ckpt 
evaluated at the scale fx B and computing the integral of <fi. By adding the symmetric term 
with fiA and /j,b interchanged, we write the contribution o\ from the terms of the first order 
in Sa e s as 
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where the function 6 is defined as 
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(12) 
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Next we use the fact that the requirement of consistency with the operator product 
expansion constrains the structure of the effective coupling variation 5a c s HH. Let us 
define the moments 



A^} = — 
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In order that the ultraviolet behavior of the running coupling be not ruined by the vari- 
ation 5a c g, the moments with n = and p integer have to vanish This implies 
that only terms in that are nonanalytic in fj,\ for small fx\ can contribute to power cor- 
rections. To find these terms, we examine the dominant region of integration in Eq. ([?]), 
fi A ^ k 2 ^ Qa Qb- By evaluating the contribution from this region, we obtain the following 
approximate expression for the function (f>: 
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where 



U 2 (z) 



z dt 
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ln(l - t) 
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Li 2 (t) 



(15) 



This formula is accurate up to the leading power in r) = ^\/{QaQb) and up to the single 
logarithms. Higher powers in rj as well as constants associated with the leading power are 
changed by the terms that we have dropped in performing the integral (|7|). To determine 
the nonanalytic behavior that controls the leading power correction to a we expand around 
r] = 0: 



<P(Q 2 A ,Q 2 B ;^ A ) 



+ 
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Thus, from Eqs. ( |TTD and ( [TB] ) we find 
(8 a £e 2 ) 2 2a PT {Q A Q B 



1 J, 3 1 1 , 2 1 
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{analytic terms} > + 0(r/ 2 ) 



C(2)~1.64. (16) 



where 
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(18) 



Here we have written oiw{QaQb\/^) = &pt(QaQb) + C(«p T \nrj) and neglected the higher 
order term on the right-hand side. In the definition of L in Eq. ([18]) we have implicitly taken 
the scale in the logarithm to be 1 GeV 2 . 
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Eq. (|T7|) gives the result for the leading power correction to the 7* 7* cross section in 

(i) 

terms of the dimensionful nonperturbative parameters A 2 . These parameters are thought 
of as being universal and should be determined by fits to experimental data. The coefficients 
of these parameters are given in Eq. (|17]) up to terms that vanish like Q 2 / s at large s. We 
see that the contributions in A*ip and are enhanced by double and single logarithms of 
the hard scale Q A Qb- 

So far, there have been attempts to study the nonperturbative moments A® based on 
data for deeply inelastic lepton-nucleon scattering and for hadronic final states in e + e~ 

annihilation. Data on the structure function F 2 suggest that A% ~ — 0.15GeV 2 fllQf . Very 

(2) (3) (2) 

little is known about A 2 and A 2 at present. A 2 enters in the power correction to the mean 

value of the three-jet resolution in e + e~ annihilation. A recent analysis of data recorded 

at PETRA suggests that this correction should be very small Assuming that the 

contribution of A^ predominates in Eq. (|TT|), owing to its enhanced coefficient, we expect 

the leading power correction to be positive. 

To complete our analysis, we need to show that, as anticipated, the term of the second 

order in Sa e s in Eq. (H does not contribute to the leading power correction. We write this 

term as 



a 2 = — / — 2^5a e fr(^i) / — T~ ^ a es(^B) • ( 19 ) 
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This contribution probes the function <fi in the corner of the phase space where both ii\ 
and fx% are small. We may restrict the k -integration in Eq. (^) to the region k ^Q 2 , with 
Q 2 being of the order of the hard scales Q A , Q 2 B , and we may rewrite the function as a 
sum of two pieces, each proportional to a logarithmic derivative with respect to one of the 
dispersion variables: 



Ai Qb'i Ai A) 
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g{k 2 /Q 2 A )g{k 2 /Q 2 B ) 
g(k 2 /Q 2 A )g(k 2 /Q 2 B : 



(20) 



The pole at fi 2 A = fi 2 B cancels in the sum of the two pieces. This may be checked, for 
instance, by substituting the expression for g, expanding it for small k 2 and performing the 
integrals explicitly. Taking this cancellation into account, the first piece gives rise to terms 
that are analytic in /j, B and possibly have nonanalytic contributions in fi 2 A . Analogously, 
the second piece gives rise to terms that are analytic in [j? A and possibly have nonanalytic 
contributions in fi B . Since, as already noted, integer moments of Sa e s vanish, for each term 
we get a vanishing contribution in Eq. ( |l9|) from either the integral in n 2 A or the integral in 
fi B . Therefore, a 2 does not contribute to the leading power correction. 

Finally, let us consider the scattering of longitudinally polarized photons. The function 
g( L ) for this case has a structure analogous to Eq. (|9p but with a different splitting function, 
p( L ) |2[] . The longitudinal splitting function vanishes at the endpoints: 

P {L) = y/2x(l -x) . (21) 
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This can be seen as being associated with the absence of aligned-jet terms for longitudinal 
photon scattering at high energies As a result, by performing a calculation analogous 
to the one described above for the transverse case, we find that the small-?] expansion of 
0( L ) has at most single logarithms: 



(Ql,Ql;v)- 8 



q / 
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(22) 



Therefore, for the power correction to the longitudinal 7* 7* cross section en L ) we get 



4 L) 



(8 a Eelf 2a PT {Q A Q B ) 
VQWb 



(23) 



We observe that the structure of the power correction is considerably more complicated in 
the transverse case. In the longitudinal case, the power correction only depends on one 
nonperturbative moment, A2 . Moreover, in the longitudinal case the coefficient is not 
enhanced by logarithms of Qa Qb- For both transverse and longitudinal scattering, it would 
be interesting to use experimental data for studying the moments and for testing the 
dispersive structure of the power-behaved terms. 
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